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Abstract 

We study the stochastic control problem of maximizing expected utility from termi- 
nal wealth under a non-bankruptcy constraint. The wealth process is subject to shocks 
produced by a general marked point process. The problem of the agent is to derive the 
optimal insurance strategy which allows "lowering" the level of the shocks. This opti- 
mization problem is related to a suitable dual stochastic control problem in which the 
delicate boundary constraints disappear. We characterize the dual value function as the 
unique viscosity solution of the corresponding a Hamilton Jacobi Bellman Variational 
Inequality (HJBVI in short). 
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1 Introduction 



We study the optimal insurance demand problem of an agent whose wealth is subject to 
shocks produced by some marked point process. Such a problem was formulated by Bryis 
[3j in continuous-time with Poisson shocks. Gollier p3] studied a similar problem where 
shocks are not proportional to wealth. An explicit solution to the problem is provided 
by Bryis by writing the associated Hamilton-Jacobi-Bellman (HJB in short) equation. In 
Bryis [3 J and Gollier |14j , they modeled the insurance premium by an affine function of the 
insurance strategy 9 = (Ot)t£[o,T] which is the rate of insurance decided to be covered by 
the agent. If the agent is subject to some accident at time t which costs an amount Z, then 
he will pay OtZ and the insurance company reimburses the amount (1 — 6t)Z. They didn't 
assume any constraint on the insurance strategy which is not realistic. 
In risk theory, Hipp and Plum [9 J analysed the trading strategy, in risky assets, which 
is optimal with respect to the criterion of minimizing the ruin probability. They derived 
the HJB equation related to this problem and proved the existence of a solution and a 
verification theorem. When the claims are exponentially distributed, the ruin probability 
decreases exponentially and the optimal amount invested in risky assets converges to a con- 
stant independent of the reserve level. Hipp and Schmidli [10J have obtained the asymptotic 
behaviour of the ruin probability under the optimal investment strategy in the small claim 
case. Schmidli [22J studied the optimal proportional reinsurance policy which minimizes 
the ruin probability in infinite horizon. He derived the associated HJB equation, proved 
the existence of a solution and a verification theorem in the diffusion case. He proved that 
the ruin probability decreases exponentially whereas the optimal proportion to insure is 
constant. Moreover, he gave some conjecture in the Cramer-Lundberg case. Hojgaard and 
Taksar [11] studied another problem of proportional reinsurance. They considered the issue 
of reinsurance optimal fraction, that maximizes the return function. They modelled the 
reserve process as a diffusion process. 

Touzi [23] studied the problem of maximizing the expected utility from terminal wealth 
when the insurance strategy is valued in [0, 1] at each time . He modeled the wealth pro- 
cess by a Doleans-Dade exponential process. He assumed a boundedness assumption on the 
jump term which guarantees the positivity of the wealth process. He solved this stochastic 
control problem by using duality method. 

Duality method was introduced by Karatzas et al. [15] and Cox and Huang [5]. Cox 
and Huang characterized the optimal consumption- portfolio policies when there exist non- 
negativity constraints on consumption and on final wealth. They gave a verification theorem 
which involves a linear partial differential equation unlike the nonlinear Bellman equation. 
In few cases they constructed the optimal control. Extensions to the case of constrained 
investment are considered by Cvitanic and Karatzas [6] and to the case of incomplete mar- 
kets by Karatzas et al. [16]. Typically, in incomplete markets, we have to solve a dual 
problem which leads in the Markov case to nonlinear partial differential equation. 
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In this paper, we model the claims by using a compound Poisson process. The insur- 
ance trading strategy is constrained to remain in [0,1]. We impose a constraint of non- 
bankruptcy on the wealth process Xt of the agent for all t. The objective of the agent is to 
maximize the expected utility of the terminal wealth over all admissible strategies and to 
determine the optimal policy of insurance. 

In our case the wealth process positivity constraint is a real one unlike the problem formu- 
lated in Touzi [23] . 

Our stochastic optimization problem is a particular case of a general structure of problems 
developed in Mnif and Pham [19] who considered the following optimization problem: 



where X + [x) := {x + X : X G X s.t X t > for all < t < T} where X is a family of 
semi-martingales. Existence and uniqueness of solution of problem (jl.ip is then proved. 
The optimal solution characterization is obtained from a dual formulation under minimal 
assumptions on the objective function. 

In this paper, we study the dual value function by a PDE approach. The dual problem 
appears as a mixed control/singular optimization problem with dynamics (Yj = ZtDt, t G 
[0,T]) governed by a classical control term Z which comes from the insurance strategy and 
a singular term D which comes from the state constraint. 

The originality of this paper is to study a stochastic control problem with state constraint. 
The wealth of the investor must be non-negative even after a jump. The duality method 
is not another alternative to solve this problem. In fact the primal problem leads to a 
HJB equation with boundary conditions. Because of the state space constraints these 
boundary conditions are not obvious to obtain. However, these delicate boundary conditions 
disappear in the dual problem. The regularity of the dual value function is not obvious 
to obtain. This explains the use of the notion of discontinuous viscosity solutions. The 
comparison theorem is not proved in a general framework since the operator which appears 
in the HJBVI contains an inf on a unbounded set which makes it discontinuous. In this 
paper, we prove the comparison theorem only when the space of claims is a finite one. 
The paper is organized as follows. Section 2 describes the model. In Section 3, we formulate 
the dual optimization problem and we derive the associated HJBVI for the value function. 
In Section 4, we prove that the dual value function is a viscosity solution of our HJBVI. In 
Section 5, we prove a comparison theorem. 

2 Problem formulation 

Let (f2, T, P) be a complete probability space. We assume that the claims are generated by 
a compound Poisson process. More precisely, we consider an integer-valued random measure 
fi(dt, dz) with compensator ir(dz)dt. We assume that ir(dz) = gG(dz) where G(dz) is a 
probability distribution on the bounded set C C JR + and g is a positive constant. In this 



max E[U(X T )], x G M, 
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case, the integral, with respect to the random measure /i(dt,dz), is simply a compound 
Poisson process: we have f f c zfj,(du,dz) = Yli^i^i, where N = {N t ,t > 0} is a Poisson 
process with intensity q and {Zi,i £ JN} is a sequence of random variables with common 
distribution G which represent the claim sizes. 

Let T > be a finite time horizon. We denote by IF = {J-t)o<t<T the filtration generated 
by the random measure fi(dt,dz). 

By definition of the intensity ir(dz)dt, the compensated jump process: 

fl(dt,dz) := fi(dt,dz) — ir(dz)dt 

is such that {/u([0,t] x B),0 <t< T} is a (P,JF) martingale for all B £ C, where C is the 
Borel ex-field on C. 

An insurance strategy is a predictable process 9 = (9t)o<t<T which represents the rate 
of insurance covered by the agent. We assume that the insurance premium is an affine 
function of the insurance strategy. Given an initial wealth x > at time t and an insurance 
strategy 9, the wealth process of the agent at time s £ [t,T] is then given by : 



■= x + {a - (3{l - 9 U )) du - / / 6 u zn(du,dz). (2.1) 
Jt Jt Jc 

We assume that a > (3 > which means that the premium rate received by the agent is 
lower then the premium rate paid to the insurer. In the literature, this problem is known 
as a proportional reinsurance one. The agent is an insurer who has to pay a premium to 
the reinsurer. We impose that the insurance strategy satisfies: 

9 S £ [0, 1] a.s. for all t < s < T. (2.2) 

We also impose the following non-bankruptcy constraint on the wealth process: 

Xl' x ' 9 > a.s. for all * < s < T. (2.3) 

Given an initial wealth x > at time t, an admissible policy 9 is a predictable stochastic 
process (9 s )t< s <T , such that conditions (|2.2|) and (|2.3p are satisfied. We denote by A(t, x) 
the set of all admissible policies and S(t,x) := {X f ' x ' such that 9 £ A(t,x)}. 
Our agent has preferences modeled by a utility function U. 

We assume that the agent's utility is described by a CRRA utility function i.e. U(x) = ^j-, 
where r\ £ (0, 1). 

We denote by / the inverse of U' and we introduce the conjugate function of U defined by 

U(y) := sup{U(x) - xy}, y>0 

= U{I{y))-yI{y). (2.4) 

A straightforward calculus shows that U(y) = where 7 = r^— and U'iy) = —I(y) for 

all y > 0. 
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The objective of the agent is to find the value function which is defined as 

v(t,x):= sup E(U(X^ X ' 6 )). (2.5) 

eeA(t,x) 

Notations: The constants which appear in the paper are generic and could change from 
line to line. 



3 Dual optimization problem 

First we introduce some notations. Let x > and t E [0, T]. We denote by ViS^tjx)) the 
set of all probability measures Q ~ P with the following property: there exists A £ I p , set 
of non-decreasing predictable processes with Aq = 0, such that : 

X — A is a Q — local super-martingale for any X E S(t, x). (3-1) 

The upper variation process of S{t,x) under Q £ P(S(t,x)) is the element A S ^ X \Q) ml p 
satisfying (Jin} and such that A - A s{t > x \Q) E l p for any A E l p satisfying (1331) . 
From Lemma 2.1 of Follmer and Kramkov |13|, we can derive V(S(t,x)) and A S ^ X \Q). 
This result states that Q € V(S(t, x)) iff there is an upper bound for all the predictable 
processes arising in the Doob-Meyer decomposition of the special semi-martingale V £ 
S(t,x) under Q. In this case, the upper variation process is equal to this upper bound. 
It is well-known from the martingale representation theorem for random measures (see e.g. 
Bremaud |2j) that all probability measures Q ~ P have a density process in the form : 

ZP = s(J S J^p u (z)-l)fl(du,dz)] , se[t,T], (3.2) 

where p £ Ut = {(p s (z)) t < s <T predictable process : p s (z) > 0, a.s., t < s < T,z £ C, 
It" So (U°SPs(^)| + p s (z)ir(dzj\ds < oo and E[Z%] = 1}. 

By Girsanov's theorem, the predictable compensator of an element X 9 £ S(t,x) under P p 



Af = I (a-/3)du+ / 9 U (J3- / p u (z) zn{dz))du. 



Z P T .P is : 



a 

We deduce from Lemma 2.1 of Follmer and Kramkov [13] that V(S(t, x)) = {P p : p £ Ut} 
and the upper variation process of P p is : 

A s s ^){pp) = (\ a -p)du+ [ Pu (z)z7r(dz))+du. 

Jt Jt Jc 

From the non-decreasing property of U, we have 

v(t,x)= sup E[U(H)], 

Hec+(t,x) 
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where C + (t,x) = {H G V\(F T ) ■ Xp° > Ha.s. for 9 G A(t,x)}. It is easy to check the 
convexity property of the family C+(t, x). For the closure property of this family in the 
semi-martingale topology, we refer to Pham [20J. 

The semi-martingale topology is associated to the Emery distance between two semi- 
martingales X 1 and X 2 defined as : 



D E {X l ,X 2 ) := ^2- n E 



n>l 



sup \Xl-Xf\Al 

0<t<T/\n 



We refer to Memin |17j for details on the semi-martingale topology. Since C+(i, x) is convex 
and closed and using the optional decomposition under constraints of Follmer and Kramkov 
|13j . Mnif and Pham [19] gave the following dual characterization of the set C+(t,x) 

HeC + (t,x) (3.3) 



J(H) := sup E 

ZeV°(t,x) ,reTt 



Z T Hl T=T - / Z u (a-/3+((3- / p u (z) zit(dz)) + )du 
Jt Jc 



< x, 



where V°(t, x) is the subset of elements P p G V(S(t,x)) such that A^S t,x \P p ) is bounded 
and l~t is the set of all stopping times valued in [0, T\. 

As a corollary (see their corollary 4.1 ), they deduce that the set of admissible insurance 
strategies .4.(0, x) is non empty iff x > b := (/3 — a)T. 

Remark 3.1 // the agent initial wealth is equal to x = b and since A(0, x) is not empty, 
then the wealth process is given by X^ ,x ' 9 = (/3 — a)(T — t) for all < t < T and so the 
only admissible strategy is 0t = for all < t < T which implies that the dynamic version 
of the value function satisfies v (i, (/3 — a)(T — t)) = for all < t < T. These boundary 
conditions obtained from the duality approach are not obvious from the primal approach. 

Now, we fix some initial wealth x > b. We make the following assumption 

Assumption 3.1 We assume that there exist 7 G (0, 1), Q G V°(t,x) with density Zt = 
^ satisfying 

(Zt)- 1 g lf(P) 

for some p > jh?- 

Under Assumption I3.lt Existence and uniqueness of problem ()2.5[) are proved (see Mnif 
and Pham pj|]). We focus now on the study of the dual formulation. 

The two following lemmas allow us to give an expression of the dual value function. We 
denote by V t the set of nonnegative, nonincreasing predictable and cadlag processes D = 
(D s )t< s <T with D t = 1, 

y°(t) ■= {Y p ' D = Z p D, Z p G V°(t,x), D G V t } 
for all t G [0,T] and L\{F T ) is the set of nonnegative J^-measurable random variables. 
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Remark 3.2 We omit the dependence of y°(t) in the initial wealth x, since x is fixed in 
all the paper. 

Remark 3.3 The set T>t is introduced in the paper of Elkaroui and Jeanblanc [8] in the 
continuous case. In our case, we must enlarge this set to include cadlag processes and 
extend some results of Mnif and Pham [19] . 

Lemma 3.1 For all x > b, X G S(t,x), Y = ZD, Z G V°(t,x), D G V t , the processes : 

Z.X - / Z u (a - f3 + ((3 - / p u (z) zn(dz)) + )du, 
Jt Jc 

and FX — J Y u (a — f3 + (/3 — f p u (z) z ir(dz)) + )du 
are supermartingales under P. 



c 



Proof. By definition of V°(t, x), the process Z.(X, — f t (a— /3+(/3— f c p u {z) z Ti(dz)) + )du) is 
a P-local supermartingale. From Theorem VII. 35 in Dellacherie and Meyer [7], the process 

M = ZX.- [ Z u {a-l3 + (/3- [ Pu (z) z 7t(dz)) + )du (3.4) 
Jt Jc 

is a P-local supermartingale. Moreover, M is bounded from below by the random variable 

— j t Z u (a — j3 + {p — f c p u (z) z ir(dz)) + )du), which is integrable under P. We deduce by 

Fatou's lemma that M is a P-supermartingale. On the other hand, by Ito's product rule 

and since D is predictable with finite variation, we get : 

Y T X T - [ T Y u (a - (3 + (P - [ p u (z)zn(dz)) + )du 
Jt Jc 

= Y t X t + j D u -dM u + ! Z u -X u -dD u + ]T AM U AD U . (3.5) 

From Equation (|3.4p . we have 

AM, / iff s = Ti = ini{u > t s.t.N u = i}, i G IN 

where Ti is the time of the claim number i and N is the Poisson process representing the 
number of claims. If the stopping time Tj is accessible, then from Theorem 15.1 in Rogers 
and Williams [21], we have Z T - = Z^ which is false and so (Tj)j g jv is a sequence of totally 

i 

inaccessible stopping times. On the other hand, D is predictable process and so from 
Lemma 27.3 in Rogers and Williams [21] we have AD T = for every totally inaccessible 
stopping time r. Since AM S = if s ^ Ti and AD S = if s is a totally inaccessible 
stopping time, we have 

AM S AD S = ds®dP a.s. (3.6) 

Since D is nonnegative and nonincr easing, and Z, X are nonnegative, this shows that the 
process : YX — f t Y u (a — j3 + (/3 — j c p u (z) zi:{dz)) + )du is a P-local supermartingale, 
bounded from below by an L 1 (P) random variable, and hence a P-supermartingale. □ 
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Lemma 3.2 For all H S L° + {Ft), we have : 

fT 



J(H) = sup E 

Yey°(t) 



Y T H 



[ Y u (a-p + (p- [ p u (z)zn(dz))+)di 
Jt Jc 



< x (3.7) 



Proof. Fix some H £ L° + {Ft)- Given an arbitrary r S 7~t, we define a sequence {D n ) n of 
elements in T>t by : 



exp 



n\ T < u du I , t < s < T, n € IN. 



We then have for all Z £ V u {t,x) : 
E ' 



Z T D T H - [ Z u Dl(a -P + (P- [ Pu(z) z n(dz))+)d 
Jt Jc 

Z T D T H - [ Z u D u (a - p + (fi - [ Pu {z) z n(dz)) + )du 
Jo Jc 

Since D™ —> l u <r a.s., for all t < u < T, we have by Fatou's lemma : 



< sup E 

zeP°(t,x),Dev t 



E 



Z T Hl T=T - [ Z u {a-P + {P- [ p u (z)zn(dz)) + )du 
Jo Jc 

Z T D T H- [ T Z u D u (a - P + (P - [ p u ( 
Jt Jc 



< sup E 

zev°(t,x),D<=v t 



z) z ir(dz)) + )du 



Identifying a probability measure Q G V°(t,x) with its density process Z, we then obtain 
from Bayes formula : 



J{H) = sup E Q 

QeV (t,x),r£Tt 

< sup E 

ZEV°(t,x),DeV t 



H1 T=T — / (a - P + (P — p u (z) zir(dz)) + )du 



C 



Z T D T H - [ Z u D u {a -P + (J3- f Pu(z) z ^{dz)) + )du 
Jo Jc 



Conversely, by the supermartingale property oiYX —J Q Y u (a— /3+(/3— J c p u (z) z Tr(dz)) + )du 
for any Y £ y°(t), see Lemma 13. 11 we have : 

fT 



sup E 

Yey°(t) 



Z T D T H- [ Z u D u {a-P + {P- [ Pu {z)zTT{dz)) + )du 
Jt Jc 



< x, \/H E C+(x)(3.8) 



Now, by characterization 
inequality (|3.8p that 



any H £ LPAFt) lies in C+(J(H)). We then deduce from 



sup E 

Yey°{t) 



Y T H- f Y u (a-P + (P- [ p u (z)zir{dz)) + )du] < J(H), 
Jo Jc 



which proves the required equality (|3.7p . 
The dual problem of (|2.5p is written as: 

fT 



□ 



v(t, v) := inf E 



U{yY£ D )+ I yY P' D (a-P + (P- I Pu (z)z7r(dz)) + )d 



c 



(3.9) 



We shall adopt a dynamic programming principle approach to study the dual value function 
(|3.9p . We recall the dynamic programming principle for our stochastic control problem: for 
any stopping time 0<T<T,0<t<T and < h < T — t, 

^ = Y J% 0{t) E H^ t + h ^ T ^<^) (3-10) 

(t+h)AT Y£ D (a - p + (p - Pu {z) z ir(dz)^J ^ du 

where a A b = min(o, b) ( see e.g. Fleming and Soner |12|). 

Lemma 3.3 Let t G [0,T] and Y P > D G y°(t). Then the process Y P,D evolves according to 
the following stochastic differential equation 



dYf> D 

with 



Y?: D ( K -dL s + f (p s (z)-l)fi{ds,dzfj, (3.11) 



dL s = -^il {Ds>0} , t < s < T, L t - = 0. (3.12) 

Proof. By Ito's product rule, we have 

dYf' D = Z s -dD s + D s -dZ s + AZ S AD S . 

From Equation (|3,2p . we have 

AZ S ^ iff s = Ti = inf{u > t s.t.N u = i}, i G IV. 

Repeating the same argument as in equation (|3.4p . we have AZ S AD S = ds g) dP a.s. 
and so equation (|3.11|) is proved. □ 

We denote by Ct the set of adapted processes (L s )t< s <T with possible jump at time s = t 
and satisfying equation (|3.12p . 

The Hamilton Jacobi Bellman Variational Inequality arising from the dynamic program- 
ming principle (|3.10p is written as 

( dv ( dv \ dv 1 

minj — (t,y) + H Uy,v,— j ,-— (t,y)\ = 0, (t,y) G [0,T) x (0,oo), (3.13) 

with terminal condition 

v(T,y) = U(y),y G (0,oo), (3.14) 

where 

H (t, y, 5, ^) := mf | A" A, y, 5, |j) + y (a - /3 + (/3 - jT p ( 2 ) * | , 
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A" (t, y, v, ^\ := (v(t, p(z)y) - v(t, y) - (p(z) - l)y^(t, 2/)) ir(dz), 

and £ := j/j positive Borel function denned on C s.t. J c ^| log p(z)\ + p(;s:)j7r(<iz) < oo|. 
This divides the time-space solvency region [0, T) x (0, oo) into a no-jump region 

( dv ( dv 

Ri = Ut,y) € [0,T] x (0,oo), s.t. -^(t,y) + H \ t,y,v, — 

and a jump region 

^2 = | (t, y) G [0, T] x (0, oo), s.t. j-(t, y) = 

4 Viscosity solution 

In this section, we provide a rigorous characterization of the dual value function v as a 
viscosity solution to the HJBVI f|3. 13j) . The function v is not known to be continuous and 
so we shall work with the notion of discontinuous viscosity solutions. We prove that the 
dual value function lies in the set of functions -D 7 ([0,T] x (0,oo)) defined as follows: 

D 7 (\0,T] x (0,oo)) :=\f: [0,21 x (0,oo) -> M such that 



(M x (0,oo)) :={/: [0, T] x (0,oo) 



\f(t,y)\ . , 
sup < oo and sup r — -r ttttt < oo > . 

y>oy + y- J x>o,y>o \x - y\{l + x-^i+ l ) + / 

The following lemma gives some properties of the dual value function v. 

Lemma 4.1 We assume that there exists a solution the the dual problem (|3.9p . The fol- 
lowing properties hold: 

1) The dual value function v is convex in y, 

2) The dual value function v satisfies the following growth condition 

su P H!l^ <00 (4.1) 

y >o y + y-~< 

3) The dual value function v satisfies 

\v(t,yx) - v(t,y 2 )\ ( s 

sup _, . —. , < oo. (4.2) 

yi>0,J/ 2 >0 \y x - y 2 \(\ + yi ^ +L > + y 2 ^ +L >) 

Proof. See Appendix. □ 



Remark 4.1 When we assume that there exists a solution the the dual problem (|3.9p , one 
can use the conjugate duality relation proved in theorem 5.1 of Mnif and Pham U9\j . 
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Since the dual value function v is locally bounded, the upper and the lower semi-continuous 
envelope of the function v are well-defined. The definition of the upper and the lower semi- 
continuous envelope of a function cp are given as follows. 

Definition 4.1 (i) The upper semi- continuous envelope of a function (p is 

<P*(t,y) = limsup 4>{t',y'), for all (t, y) G [0, T) x (0, oo). (4.3) 

(V ,»')->(*,») 

t'6[0,T),y'>0 

(ii) The lower semi- continuous envelope of a function (p is 

Mt,y)= liminf <P{t' ,y'), for all (t,y) G [0,T) x (0, oo). (4.4) 

t'6[0,T),y'>Q 

Since the continuity of the Hamiltonian in his arguments is not obvious, we define lower 
semi-continuous envelope of H by 

fl»(t,y,^,&:= liminf F(t', y', ^, ^) for all (i, y) G [0, T) x (0, oo). 

oy {t> , y >)^(t,y) ay 

t'e[0,T),y'>0 

Adapting the notion of viscosity solutions introduced by Crandall and Lions [1] and then by 
Soner [23jfor first integrodifferential operators, we define the viscosity solution as follows: 

Definition 4.2 (i) A function cp is a viscosity supersolution of \3.13\) in [0, T) x (0, oo) if 
min(^(t,y) + #* ft, y, ^ ^ ) , -^-{t, y) ) < 0, (4.5) 



at \ oy J oy 

whenever ip G C 1 ([0, T] x (0, oo)) and cp* — ip has a strict global minimum at (i,y) G 
[0,T) x (0,oo). 

(ii) A function cp is a viscosity subsolution of f3.13\) in [0, T) x (0, oo) if 

m i„{f(f, S) + H(f, S ,*,|),-|(f, 5 )}>0, (4.6) 

whenever ip G C 1 ([0,T] x (0,oo)) and cp* — ip has a strict global maximum at (i,y) G 
[0,T) x (0,oo). 

(Hi) A function <p> is a viscosity solution of h3.13\) in [0, T) x (0, oo) if it is both super- solution 
and sub-solution in [0,T) x (0, oo). 

Remark 4.2 As it is seen in the definition of viscosity solutions, we use the lower semi- 
continuous envelope of the Hamiltonian H. In fact to prove that the dual value function v 
is viscosity solution of our HJBVI, we need only the regularity of to derive inequality 

(Mp . 

The following theorem relates the dual value function v to the HJBVI (|3.13p . 

Theorem 4.1 The dual value function v is a viscosity solution of h3.13\) in [0, T) x (0, oo). 
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Proof. See Appendix. □ 

The HJBVI (|3,13|) associated to our problem does not provide a complete characteriza- 
tion of the dual value function v. We need to specify the terminal condition. From the 
definition of v, it's obvious that v(T,y) = U{y). Since we use the notion of discontinuous 
viscosity solutions, we need to characterize v*(T,y) and v*(T,y) for all y G (0,oo) which is 
the object of the following theorem. 

Theorem 4.2 The terminal conditions of the upper and lower semi- continuous envelope 
of v satisfy the following inequalities 

v*(T,y) := limsup v{t',y') < U(y), for all y G (0,oo), (4.7) 

t'£[0,T),y'>0 

v*{T,y) : = limsup v(t',y') > U(y), for all y G (0,oo). (4.8) 

(t' ,y')^(T,y) 

t'e[0,T),y'>0 

Proof. We first prove inequality (|4.7p . Suppose on the contrary that there exists a constant 
rj > such that v*(T,y) > U{y) + 2rj. From the definition of v* , there exists a sequence 
((tn,y n ))neiN such that (t n ,y n ) — > (T,y) and v(t n ,y n ) — > v*(T,y) when n tends to 
infinity, which implies 

r-T 

rpn,D n \ I Vpn,Dr, 



E[U{y n Y^) + / y n Yf^(a - p + (/3 - / Pns {z) ZTr(dz))+)ds] > U(y) + V , 

for all YP^ Dn G y°(t n ) .Choosing p ns = 1 and D ns = 1 for all s G [t n ,T], we obtain 

U(y n ) + I y n {a-l3 + {!3- ! z Ti{dz)) + )ds > U(y) + r). 

Sending n to infinity, we have U{y) > U(y) + rj which is wrong and so inequality (|4.7p is 
proved. 

We prove now inequality (|4.8p . From the definition of S#, there exists a sequence ((i n , y n ))n£]N 
such that (t n ,y n ) — > (T,y) and v(t n ,y n ) — > C*(T, y) when n tends to infinity. 
From the definition of the dynamic version of the value function, we have for all x > 
(P — a)(T — t n ) 

v(t n ,x)>U(x + (a-(3)(T-t n )). (4.9) 

Let e > 0, there exists uq G IN such that for all n > no, we have 

U(x + (a-P)(T-t n )) >U(x)-e. (4.10) 

Using the conjugate duality relation of Theorem 5.1 of Mnif and Pham [19] and relations 
3]) and dH0|), we obtain 

v(t n ,y n ) = max [v (t n , x) - xy n ] 

x>(/3-a)(T-t n ) 

> max \U(x) — xy n ] — e. (4.11) 
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Since (t n ,y n ) — > (T,y) when n tends to infinity, there exists n\ E IN such that for all 

n > ni, I(y n ) > ((3 — a)(T — t n ) and so max W( x )~ x Vn] = U(y n ). For n > uq Vrti, 

x>(P-oi)(T—t„) 

inequality (|4.1ip implies 

v{t n ,y n ) > U(y n ) ~ e. 

Sending n to infinity and e to 0, we prove inequalitv ([4.8p . □ 



5 Uniqueness 

We turn now to uniqueness questions. Our next main result is a comparison principle for 
discontinuous viscosity solutions to the HJBVI (|3. 13[) . It states that we can compare a 
viscosity sub-solution and a viscosity super-solution to the HJBVI ()3.13|) on [0, T) x (0, oo), 
provided that we can compare them at terminal date as usual in parabolic problems. 
The main difficulty in the comparison theorem comes from the discontinuity of the Hamil- 
tonian. Here, we assume that the claims take values in the set C = {6\, 62, S^}, <5j > 0, 
1 < i < d. In this case the Hamiltonian contains an inf on a bounded set which makes it 
continuous. We denote by m, 1 < i < d the intensity of the Poisson process associated to 
the claim having the size <5j. The set £ is defined as follows : 

£ := {p = (pi)i<i< d , Pi > 0, 1 < i < d} . (5.1) 
The Hamiltonian H is given by 

where 

The following lemma sates some properties of the functions t>* and v*. 

Lemma 5.1 We assume that there exists a solution the the dual problem (|3.9|) . The fol- 
lowing properties hold : 

1) The functions v* and are convex in y. 

2) The functions v* and v* are nonincreasing on (0, 00). 

3) The functions v* and v* satisfy the following growth condition 

\v*(t,y)\ \v*(t,y)\ 
sup , sup < 00 (5.3) 

y >o y + y 7 y >o y + y 7 

4 ) The functions v* and v* satisfy 



\v*(t,yi)-v*(t,y 2 )\ 
sup —, . —. . < 00. 

yi>0,j/ 2 >0 \ yi - 2/aKl + y x i7+ij + y 2 ^ 7+ij ) 
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\v*(t, yi )-v*(t,y 2 )\ 
sup —. , —. , < oo. 

yi>o m >o \yi-y 2 \(l + y 1 17 j + y 2 17 j ) 

Proof. 1) Fix A G (0,1) and (y,y',y") G (0,oo) 3 such that y = Ay' + (1 — A)y". From 
the definition of v*(t,y), there exists a sequence ((t n ,yn)) n such that v(t n ,y n ) — )• v*(t,y) 
when n goes to infinity. We set y n = Ay' + (1 — A)y n , then we have y' n — > y when n goes 
to infinity. Since v is convex in y, then we have 

v{t n ,y n ) < Xv(t n ,y) + (1 - A)£(i„,y^) (5.4) 
Sending n to infinity, inequality (|5.4p implies 

v*(t, y) < Xv*(t, y) + (1 - X)v*(t, y"), (5.5) 
which is the desired result. 

We turn to the convexity of v*. From Theorem 5.1 of Mnif and Pham ID] . we have 

v(t,y) = max [v(t,x)—xy]. (5-6) 

x>(f3—a){T—t) 

We fix (t,y) G [0, T] x (0, oo). From the definition of 5*, there exists a sequence ( (t n , y n )) n 
such that v(t n ,y n ) — > v*(t, y) and (t n , y n ) — > (t, y) when n goes to infinity. From equation 
()5.6|) . we have v(t n ,y n ) > v(t n , x) — xy n > v*(t, x) — xy n for all x > (/3 — a)(T — t n ). Sending 
n to oo, we have v*(t, y) > f*(t, x) — xy for all x > (/3 — a)(T — t) and so 

v*(t, y) > max [vJt,x) — xy]. (5.7) 
x>(j3-a)(T-i) 



From Theorem 5.1 of Mnif and Pham [19], there exists x such that 

v(t,y) =v(t,x) - xy (5.8) 

dv 

and Tr-(t,y) = —x. We consider a sequence ((t n ,y n )) such that v(t n ,y n ) — > vJt,y) and 
ay n 

(t n ,yn) — > (t,y) when n goes to infinity. From equation (|5.8p . there exists i n such that 

v{t n ,y n ) = v(t n ,x n ) - x n y n . Since u G £> 7 ( 0,T x (0,oo)), we have x n := |^-(t„,y n )| < 

dy 

K(l +yn il+l)) ) where X is a positive constant independent of n. Since y n — > y when n 
goes to infinity, the sequence (x n ) n is bounded, and so long a subsequence x n — > x when 
n goes to infinity and so we have 

lim v(t n ,x n ) = v*(t,y) + xy = vJt, x) (5.9) 

n — >oo 

From inequality (|5.7p and equation (|5.9p . we deduce that 

v*(t,y) = max [^(t,x)— xy]. 
x>(/?-a)(T-t) 
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and so is convex in y. 

2) Fix (t,y) G [0, T) x (0, oo). Since at time t, only D £ D ( could make a jump, we have 
y = Yt > Yj+ . From the dynamic programming principle (|3.10|) we have 

v(t, y) < v(t, y(l - 5)) for all < 5 < 1, 

and so the dual value function v in non-increasing with respect to y. This yields that v* 
and v* are nonincr easing. 

3) Fix (t,y) G [0, T) x (0, oo). From the definition of v*(t,y), there exists a sequence 
{{tn,Dnj) n such that v(t n ,y n ) — )• v*(t,y) when n goes to infinity. From the growth condi- 
tion of v, we have 

\v(t n ,y n )\ < K(y n + y~' r ), 

where K is a positive constant. Sending n to infinity, we obtain the desired result. We use 
• ., , \v*(t,y)\ 

similar arguments to prove that sup < oo. 

y>o y + y 7 

4) We prove only the first inequality. The second one is obtained by using similar arguments. 
Fix (t,yi,y2) G [0,?") x (0, oo) x (0,oo). From the definition of v*(t, 2/2), there exists a 
sequence ((i n , y2,n)) such that v(t n ,y2 >n ) — > v*(t,y2) when n goes to infinity. From 
inequality f)4.2[) . we have 

v(t n ,yi) -v(t n ,y 2 , n ) < K\yi -y 2 ,n|(l + yi iri+1) +y^i 7+1) ), 
where K is a positive constant. Since liminf v(t n , y%) > t5*(t, yi), we obtain after sending n 

n — s-oo 

to infinity 

-v*(t, y 2 ) < K\ yx - y 2 \(l + yr (7+1) + y 2 (7+1) ) - 2/i)- 

Using similar arguments, we deduce the inverse inequality and so 

\v*(t,y{) - v*{t,y 2 )\ 



SUp — ; ■ —. ■ < OO. 

yi>o,y 2 >o |yj - + y x i7+ij + y 2 (7+ij ) 



□ 



Remark 5.1 One could prove the monotinicity of v* by using viscosity solutions argu- 
ments. In fact for each e > 0, we define W(t,y) = v(t,y) — ey, (t,y) G [0,T) x (0,oo). The 

dW 

function W satisfies in the viscosity sense — — < — e, i.e. for all (yo, t/j) G ((0, 00), C ([0, T] x 
(0, 00))) such that 

(W*-iP)(t,y )= max (W*-^)(t,y) (5.10) 
ye(o,oo) 

dip 

and so we have -77— (t, yo) < — e. This proves that ^(t, .) is nonincreasing on a neighborhood 
ay 

V(y ). Let (2/1,2/2) G V(2/ ), we want to prove that W*(t,y 1 ) > W*(t,y 2 ). 
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Suppose that W*(t,y{) < W*(t, 2/2)- We consider the function V(t, y) = W*(t, y{) which 
solves 

dV 

— =0 on (2/1,2/2), (5.11) 

together with the boundary conditions V(t,y\) = V(t, 2/2) = W*(t,yi). Since W is a 
viscosity subsolution of Equation (j5.11j) . From the comparison theorem ( see Barles [I], 
Theorem 2.7 in the case of continuous viscosity solutions), we have 

inf (V - W*)(t,y) = mf((V - W*){t,y x ), (V - W*)(t,y 2 )) = 0, 

mm] 

which implies V(t,y) > W*(t,y) for all y £ [2/1,2/2]- Since W*(t,y ) < V(t,y ) = W*(t,yi) 
and ip(t,y ) > ip(t,yi), then 

which contradicts f|5. 10|) . Sending e — > + , we obtain the desired result. 



Now, we are able to prove the following comparison principle : 

Theorem 5.1 Let v\ (resp V2) £ -D 7 ([0,T] x (0,oo)) be a viscosity subsolution (resp su- 
persolution) of (|3.13p in [0, T] x (0,oo) such that vl(T,y) < V2*(T,y). We assume that v 2 * 
is convex and nonincreasing in y, then 

v*(t,y) < v 2 *(t,y), for all (t,y) G [0,T] x (0,oo) (5.12) 

Proof. See Appendix. □ 



By combining the previous results, we finally obtain the following PDE characterization of 
the dual value function. 

Corollary 5.1 We assume that there exists a solution the the dual problem (|3.9p . The 
dual value function v is the unique viscosity solution of (|3. 13[) with terminal condition 
v(T,y) = U(y) in the class of functions Z? 7 ([0,T] x (0,oo)). 

Remark 5.2 In this remark, we formally discuss the numerical implications of Corollary 
15.11 We can solve numerically the associated HJBVI by using an algorithm based on policy 
iterations. Then thanks to a verification theorem, We characterize the optimal insurance 
strategy by the solution of the variational inequality. These results are the object of the 
paper Mniffltff. 
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6 Appendix 



6.1 Proof of Lemma 14.11 

1) Prom Theorem 5.1 of Mnif and Pham [19], we have v(t,y) = max [v(t,x) — xy\. 

x>(/3—a)(T—t) 

The convexity property of v in y holds since it is the upper envelope of affine functions. 

2) Since the controls p s = 1 and D s = 1, s£ [t, T] lie in lit x T>t, we have 

v(t,y)<U(y) + Ky, (6.1) 

where K is a constant. 

Let (Z n := Z pn ,D n ) be a minimizing sequence of v(t,y). From the definition of these 
minimizing sequences, there exist e n and uq £ JN such that e n — > when n — > oo and 
for all n > no, we have 



v(t,y) > E U(yZ%DZ) 



+ yE 



(' T ZZD n u {a-f3 + {P- [ p n u {z) 
Jt Jc 



z ir(dz)) + )du 



(6.2) 



Since e n — > when n — > oo, there exists n\ G IN such that for all n > n\, we have 
£n < U(y) + y- We recall That U(y) > U(0 + ) > and so U(y) + y > since y > 0. Using 
the boundedness of D n , Jensen's inequality and the martingale property of Z n , we have: 



E 



U(yZ^)\ > U(yE[ZZ}) 

> u{y)- 



(6.3) 



For the second term of the r.h.s of inequality (|6.2j) . since D™ < 1 for all s G [t, T], using 
Fubini's theorem and the martingale property of Z n , we have 



E 



f T yZ™D n u {a-P + {P- [ Pu (z)z7r(dz)) + )du 
Jt Jc 



> y(a-f3)E 



ZlD n u du 



> y{a-P) J E\Zl\du 

> K'y, 

where K' is a constant independent of y. Inequalities (j6.3j) and (j6.4j) imply that 

v(t,y) > U(y) + K'y. 

From inequalities (|6.ip and (|6.5p . we deduce that 

\v(f,y)\ 



sup ■ 

y>° y + V 



< oo 



3) Using the convexity in y of the dual value function v, we have 
\v(t,yi) - v(t,y 2 )\ , m , i~/ m 

i ; < + wife) > 

yi - 2/2 



(6.4) 



(6.5) 



(6.6) 



(6.7) 
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where v' d is the right-hand derivative with respect to the variable y. Let {Z n := Z pn ,D n ) 
be a minimizing sequence of v(t,yi). Let 5 > and (Z' n := Z' p " ,D' n ) be a minimizing 
sequence of v(t, y\ +5). From the definition of these minimizing sequences, there exist 
e f n and no G IN such that e n — > 0, e f n — > when n — > oo and for all n > no, we have 



+ yiE 



7 



2%DZ(a-p + {p- / p n u (z)z7r(dz)) + y u 



(6.8) 



and 



v{t, yi + S) > { yi + 5)^E 
+ (yi + S)E 



(Z'™D'f)-~* 



Z':D':(a-P + {P- f p':(z)zn(dz)) + )du 
J c 



e' 



(6.9) 



Using the definition of v(t,y\) and v(t,yi +5), we have 



v(t, yi ) < y^E 
and 



(ZfpD%)-i 
7 



+ yiE 



Z?D*(a -P + tf- J f%(z)z*{dz))+ )du 



v(t, yi + 5) < {yx + Sy^E 

+ (yi + S)E 



7 



Z n u D n Ja 



(6.10) 



(3 + {P- j c P n u {z)zK(dz)) + ) 



du 



Using inequalities (|6.8j) and (|6.10p . we deduce that 
v(t,yi + 5) - v{t,yi) 



< (yi + 6)^-y^ E 



+ E 



£ Z n u D n u (a-f} + (p- j c P n u {z)z Tr(cte))- 



du 



+ 



From inequality (|6.3p . we have 



E [(Z2D%)-1] > if, 



(6.11) 



(6.12) 



where if is a positive constant independent of y and 5. S ince c n — y when n — y oo, 
there exists ni such that for all n > m, we have % < 1. For the second term of the r.h.s 
of inequality (|6.1ip . since D s < 1 for all s G using Fubini and martingale property of 
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Z, we have 



E 



f Z™D n u (a-P + {f3- ( P n u (z)z7T(dz)) + )du 
t J c 









^jdu 


< aE 





< a J E\Zl\du 



< K. 



(6.13) 



Using inequalities (fU7TT|) . ([6TT2]) and (fBTTBl 

v(t, yi + 5)- v(t, yx) ( yi + 5) 



< K( 



+ D 



Sending 5 — > + , we have 



i' d (t, yi )<K(-y-^ + l). 



.14) 



Similarly, we obtain that 

v(t,yi + 5) - v(t,yi) 
5 



> (yi + 5)^-yp E 



(ZtpD!p)-"< 



7 



+ E 

We define ~p s :=C~ p 



Z?D*(a-P + (J3- [ p':(z)zn(dz)) + )du 



n 

8 ' 



(6.15) 







using the martingale property of Z p , we have 



and D s = 1 for all s G [t,T]. From the definition of v and 



v(t,yi +5) 
< { yi+ 5)-^E 

= ( yi +5)-'E 



+ (yi + S)E 



ZPD u la-f3 + (f3- / p u (z)zn(dz)) + )du 



c 



(4)" 7 

7 



+ { yi + S)(a-/3)(T-t). 



(6.16) 



From inequality (|6.9p . using the martingale of Z' n and since < D'J 1 < 1 for all s G [i, T], 
we have 



7 



t5(t,Ift + *) > (yi + <5)-^ 
from (|6.16p and (|6.17|) . we deduce that 

E [{ZpD!p)-"<} < E Uzlj,) 
We know that Z p is given by the formula 

Z£ = exp(g(t - T){C P - l))(CP) NT ~ Nt 



+ ( yi + 5){ a -f3)(T-t)-e' n (6.17) 



1 The solution of the SDE dZ t = Z t -dH t is given by the Doleans Dade exponential formula Z t = 
exp exp (—AH 3 ) where < H c > is the quadratic variation of the continuous 

0<s<t 

part of H . 
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and so E 



< oo. Since e' n — > when n — > oo, there exists ri2 such that for all 
n > ri2, we have E [(Z!p D!p)~~ ( ] < K', where K' is a positive constant independent of y\ 
and 5. Using the boundedness of D, inequality (I6.15P implies 



v(t,y 1 +6)-v(t,yi) >K ,((yi + $) 7 - V\ 7 



Sending 5 — > + , we have 

^, yi )>-^W 7+1) + i)- 

Using ([EZD, (fpTTID and (fpTTSD . we obtain 

|£(£,2/i) - ^(i, 2/2)1 



sup 



J/i>0,j/ 2 >0 | yi _ y2 | M + y~ (7+1) + ^ {7+1) 



1 . 



< OO 



.18) 



□ 



6.2 Proof of Theorem I4TT1 



We first prove that v is a viscosity sub-solution of (|3,13p in [0, T) x (0,oo). 

Let (t, y) G [0, T) x (0,oo) and ?/> G C 1 ([0, T] x (0,oo)) such that without loss of generality 

= (v* —ip)(t,y) = max (v*—ip). 

[0,T)x(0,oo) 

From the definition of v* , there exists a sequence (t n ,y n ) £ [0, T) x (0, oo) such that 

(*n>2/n) — >■ (*,2/) and v(t n ,y n ) — ► u*(i,y) when n — > oo. 

For 7] > 0, /3 ns = p a positive Borel function and _D ns = 1 for all s > t n , we set 

n := infO > t n such that (t,y n Yf n ' Dn ) $ B v (t n ,y n )} A T, 

where B v (t n ,y n ) = {(t,y) € [0,T] x (0, oo) such that \t - t n \ + \y - y n \ < r)}. By the 
right continuity of the paths, we have 9 n > t n a.s. For all < h < T — t n , the dynamic 
programming principle 



v(t n ,y n ) 



inf E 

(t n +h)A9 n 



vUt n + h) A0 n ,y n Y ( lf +h)Adn 



+ 



y n Yf> u ( a - (3 + ( f3 - p s (z) z n(dz) 



ds 



implies 



ln + i>(t n ,y n ) < E^((t n + h)A0 n ,y n Y^ )Ae ^ 

-(t n +h)A9 n / / 

+ I y n Y s pn ' Dn (a-P+[P- I P n{z) z 7r(dz) 



3.19) 



c 



ds 
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E 

+ E 
+ E 



where the sequence 7„ := v(t n ,y n ) — ip(t n , y n ) is determinist and converges to zero when n 
tends to infinity. Applying Ito's formula to ij){t n + h,y n Y t p ^ Dn ), we get 

~l r(t n +h)/\6 n 

hj t ~dt 

I r(t n +h)A6„ 

'tn 

I r{t n +h)/\6 n 



{s,y n Yf»> D ») + A* (^ s ,y n Yf"> D ",ij,^j ds 

n / ft \ \ 

j-jf VnY s Pn ' Dn \a-P+ [P- J c p n {z)zir{dz)j j ds 

Y> (s, Pn{z)y n Y p s - ' Dn ) - 4> (s, y n Y^' Dn ) p(ds, dz) 



> 



In 



By the martingale's property, the third expectation on the left hand-side of the last in- 
equality vanishes and so we obtain 



E 



1 r(t„+h)A0„ a 1 / a 1.- 

1 P ( s ,y n Yf^)+A^ ^ s ,y n Yf^,^-^ ) <1, 



h 



dt 



I r(t n +h)A8 n 
+ hi 



y n Yf"> D " (a-(3 + ({3- ^ Pn (z)zir(dz)). 



ds 



> (6.20) 



From the definition of j n , two cases are possible: 

★ Case 1: if the set {n > : 7„ = 0} is finite, then there exists a subsequence renamed 
{ln)n>o such that 7n 7^ for all n and we set h = ^/Tn- 

-k Case 2: if the set {n > : 7„ = 0} is not finite, then there exists a subsequence renamed 
{ln)n>o such that 7n = for all n and we set h = rT x . 

In 

In both cases > as n tends to 00 . We now send n to infinity. The a.s. convergence 

h 

of the random value inside the expectation is obtained by the mean value Theorem. Since 
j c n(dz) < 00 and using the definition of 9 n , the random variable 

1 f(t n +h)Ad 

l J t (s, ynYf Dn ) + A p " (s, y n Yf"> D ",i;, 

+ y n Yf Dn (a-P+tf-J^ p ns (z) z ir(dz)) + ^J ds 

is essentially bounded, uniformly in n, on the stochastic interval [t n , (t n + h) A 9 n ]. Sending 
n to infinity, it follows by the dominated convergence theorem 

^(t, y) + A?(t, y, V>, ^) + y (a - P + (0 - jf p(z) z n(dz)) + ^J > 0, 



for all p G S and so 



It remains to prove 



dib , , dib , 

-£(t,y) + H(t,y,i;-£-)>0. 



- B ^(t,y)>°. 

dy 



(6.21) 
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For 5 G (0, 1), we set 



L s 



if s = t n 
6 if s > t n , 



and so the process D is a constant for all s > t n . We choose p s {z) = p(z) for all s > t n 
and z G C, where p £ E. We have y n Y p ' D G 3^°(in)- Sending n to infinity, we have 
yY^ D = y(l — 8). Sending n — > oo in (|6.19p . by the dominated convergence theorem we 
get 

W,y)<iKt,y(l-8)). 

Sending 5 — > + , we obtain 

~^(t,y)>0. (6.22) 
dy 

Combining (|6.21j) and (I6.22h . we conclude that v is a viscosity subsolution. 

For supersolution inequality (|43|) . let ip G C x ([0,T] x (0,oo)), G [0,T) x (0,oo) such 

that (u# — , 0)(t, y) = min(t;* — ^) = 0, we need to show 

1 at + * { t,y, ^dy~) - ' (6 3) 

Suppose the contrary. Hence the left-hand side of (16.23H is positive. By smoothness of ifj 
and since H* is lower semi-continuous, there exist r/ and e satisfying: 

. j d^(t,y) ( di; \ 

mm \~d^ + * { ,y ^ , dy~) , ~dy~ {,y) j - £ ' ( } 

for all (t, y) G B v (t, y), where B v (i, y) = {{t, y) G [0, T] x (0, oo) such that \t-i\ + \y-y\ < 

7]}. By changing rj, we may assume that B v (t,y) C [0, T) x (0,oo). 

Since (i, £/) is a strict global minimizer of v* — ij) , there exists £ > such that 

min _ (v* - ip)(t,y) = £, 
(t,y)edB v (t,yy 

which implies v*(t,y) > £ + Tp(t,y) for all (t,y) G dB v (i,y) the parabolic boundary of 
B v (i,y). From the definition of v*, there exists a sequence (t n ,y n ) G [0, T) x (0, oo) such 
that y n ) — )• (t, y) and u(t n , y n ) — > y) when n — )• oo. We suppose that (t n , y n ) G 
B„(i,y). Let g 3^ (i n ) be given and the stopping time 9 n defined by 

9 n = inf{t > t n such that (t, y n Y t p ' D ) <£ B v (i, y)} A T. 

Since the control L G is singular with possible jump at t = t n , the couple (t,y n Y t p,D ) 
might jump out of B n (t,y) at t n . If the control D makes alone the latter couple jump out 
of B v (t,y), we set 9® := 9 n else 9® '■= T. In this case, the process Y decreases. We know, 
from the dynamic programming principle (|3.10p that 

v(t, y) < v(t, y(l - 5)) for all < 5 < 1, 
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and so the dual value function v in non-increasing with respect to y. However, the point 
Poisson process could contributes to the jump out of B„(i,y). In this case the dual value 
function v is not necessarily non-increasing in the direction of the jump. The control p could 
also contributes to hit the boundary of B„(t,y). To overcome this problem, we introduce 
6j the first time after t n the state process Y jumps because of the point Poisson process 
and we set On := n when y n Y p,D jumps out B v (i, y) because of the control p else On := T. 
We set also P := On A Oj. Note that, by right continuity of the paths, we have P > t n a.s. 

%■ 



Let be the stopping time defined as follows := 0® A P . 



-k On the set {0® < p }. Let (0®,y) be the intersection between dB v (i,y) the parabolic 
boundary of B v (t,y) and the line between (0® ,y n Y e D-) and (0® ,y n Y g r>) . From (|6.24p . we 
deduce that ip is non- increasing along this line in B„ (i, y) . Since the dual value function v 
is non-increasing with respect to y, we have 

v(0°,y n Y e n) > v(0%,y) > W5,v) + Z> C + £ 
Using the inequality above and applying Ito's formula to ^(t,y n Y t p ' D ), we have 

v(0°,y n Y e n) (6.25) 

> ^0°,y n Y e n-)+ti 

> 1>(t n ,y n ) + ^(s,ynYf> D )+AP(s,y n Yf> D ^,^)ds- J"" ^( s ,y n Yf' D )y n Y^ 'dL s 
+ f n (i> (s, Ps {z)y n Y^ D ) - $ (a, y n Y p l D ) fi(ds, dz) + £. 

For t n <s<0®, (16321) implies: 

%(s,y n Yf> D ) + A p ( S ,y n Yf> D ,^) 

+ y n Y s P> D (a-p + tf- J Pe (z)zn(dz))+^ > 0, (6.26) 

-?±(s,y n Yf> D )>0. (6.27) 
dy 

Substituting (16361 and (16371) into (1635]) . we have 

v{0°,y n Y e f) + y n Yf> D ( a - (3 + (f3 - [ p s (z) z n(dz))+) ds (6.28) 

> 



t n \ J c 



i>(tn,y n ) + J ip (s, p s (z)y n Y s P l D ^j - ip (s, y n Y s P l D ^j p{ds, dz) + £. 
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* On the set > 9 p }, we have 



v(9 p ,y n Y 6p ) > i>(0 p ,y n Y 6p ) (6.29) 
> ^(t n ,Vn) + J' ^(s,y n Y p > D ) + A p (s,y n Y p > D ,^, ^-)ds 

lp ^±(s,y n Y^ D )y n YPl D dL s 



+ 

tn J c 

For t n <s<6 p , ([BT24")) implies: 



/ i> (s, p s (z)y n Y s P l D ^j - ij) (s, y n Y s P l D ^) K ds > dz). 



^(s,y n Yf> D ) + A»(s,y n Yf> D ,i,,^) 



+ !hl Y a p ' D ( o - .* + ( f~ J p 8 [") :,tU/:))_ ) > <. (().:]()) 



-^ (s , yn y/^)>o. 

Substituting iO0| and (16311) into (j!T29"]) . we have 



«(0p, yn^p) + / yn^/' D ( O - P + W - I Ps{z) Z 7T(dz)) + ) ds 



c 



(6.31) 



(6.32) 



> ip(t n ,y n )+ / (i> (s,p s (z)y n Y p l D ^J - tp (s,y n Y p l D ^p(ds,dz) + e(9 p - t n ). 



Putting the two cases (|6.28p and (|6.32p together, we get 



E 



> E 



(9, y n Yt D ) + / y n Yf' D [ a - p + (/3 - / p s {z) z Tt(dz)) + ) ds 



' t n \ J c 

r 6 D / 



(6.33) 



Hes<e p} \v(9n,ynY e p f) + Jf ^F/^ (a- + (/3 - J^z) zv{dz)) + j ds 



E 



HeS>e P } Wp> Vn Y e P ; ) + J VnY p ' D \a-P + (J3-j Ps {z) z ,(dz) )+ ds 
> ^(t n ,y n ) + ZP(e% <e p ) + eE[l {g n>g p} (O p -t n )} 
Suppose that for all £' > 0, there exists Y G 3 ;0 (in) such that 

£P(0£ < 9 p ) + eE[l {e oy 6p} {e p - t n )\ < e' (6.34) 

Since # p > i n a.s. and < eE[l^ d D >Sp j(9 p — t n )] < £ , for £ sufficiently small, we deduce 
that 9% < 9 p a.s. Inequality (|6.34p implies £ < for £ sufficiently small which is false and 
so there exists £ > such that for all y G 3^ (in)? we have 



£P(#f < 9 P ) + eE[l {e n> 0p] {e p - t n )\ > (. 



(6.35) 
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Inequalities (|6.33p and (|6,35p imply 

v(t n ,y n ) + 5 n 



(6.36) 



< E 



hVnYt D ) + / y n Y s p ' D (a-/3 + (l3- / p s (z)z7c(dz)) + y s -C 
where 5 n := tp(t n ,y n ) -v(t n ,y n ). 

Since 5 n = ip(t n ,y n ) — ip(i,y) + £*(i, y) — v(t n ,y n ), there exists no £ IN such that for all 
n > no, 5 n > — |. Inequality (|6.36|) implies 

v{t n ,y n ) 

< inf E \v(9,y n Y 9 p ' D )+ f y n YP' D (a - f3 + {{5 - [ p s (z) z ir(dz)) + )ds 

which is a contradiction with the dynamic programming principle and so we conclude that 
the dual value function v is a viscosity super-solution. □ 



c 

2' 



6.3 Proof of Theorem [EH 

For e, A, S, ( > 0, we define $ : [0, T) x (0, oo) x (0, oo) — > M U {-co} as 
®(t,yi,y 2 ) := vl(t,yi) - v 2 *(t,y 2 ) - -(yi - y 2 ) 2 



8 exp (A(T - t )) ( 2/ 7+1 + % 7+1 ) - C(^r + -^T: 

v 7 y/ y 2 



Since 6J, 6 2 * E £> 7 ([0,T] x (0,oo)), there exists (^«AA,C ; X «,W ; y «,5,A,C) G [ ,T] x (0,oo) x 
(0, oo) which maximizes $. By using the inequality 

2^(t* e ' 5 ^,y* 1 e ' SAC ,y* 2 e ' 5 ' X ' C ) 

e,5,A,C , *e,<5,A,C *e,<5,A,C 



we have 



+ £ 2 *(t* £ ' w ,2/r w ) - ^(t* e ' 5 ' A ' c ,yr ,5,A ' c )- 

since v 2 * G L> 7 ([0, T] x (0, oo)), we have 



*e,<5,A,£ *e,<5,A,f 



2/2 



< C(l + 



+ 



(y«'^^)7+l (yr' 5 ' A ' C )^ +1 



-)• 



(6.37) 



Using the inequality $(i* e »W yT' S ' X ' ( , y 2 eAX ' ( ) > 1, 1) and since v\, 2 * € D 7 ([0, T] x 
(0, oo)), we have 



+ c 



1 1 



(6.38) 



+ 



^*e,5,A,C^ 7+ i (y* e ' <5 ' A ' C )*7+ 



< ^,C( 1 + ^W + ^,A,C + _J_ + 1 



(2/r ,5 ' A,c ) 7 (y 2 * e,<5,A,c ) 7 
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where C 6 ^ is a constant depending only on 5 and £. Inequality (|6.38|) implies either 

5 ( *e,6,x,(y+l + c 1 < C 5 ^ ( 1 + y*"' 5 '^ + i-^- ) (6.39) 



or 



< 5( 2/ -.^ jC)7+1 + C ~ 1 < c*C ( i + ^AA,c + 1 



(y « aa,c )7+ i - ^ ( y r ,5,A ' c )^ / ' 

Assume the first case, then there exist , > depending only on 5 and £ such that 
M 8,c < y *6,5,A,C < M 5,C_ Using i nequa lity Q8| . we obtain 

and 



which implies that < y^' 5 '^''' < . Since y* 6 ' 5 ''^ and y^ )&)X ^ are bounded from 

below, inequality (I6.37P implies 

\yT' 5M -y* 2 e ' 5M \ <C ie , (6.40) 

where C\ is a positive constant independent of e. Using the boundedness of y* e,<5,A '^ and 
y*e,5,A,£ an( ^ (|g 4Q|) along a subsequence (i* <E ' 5 ' A '£, y 1 e ' <5,A ' < ', y^'^^^) converges when e — > 0. 
Let's denote (t* 5 ' A ^, y* 5 ' A ' C , y*- 5 ^) its limit. 
From the definition of t* e ' 5 ' A >f ; two cases are possible: 

★ Case 1: If the set {e > : t* e > 5 ' A ^ = T} is not finite, then there exists a subsequence 
renamed (t* <E,<5,A '^) e such that t* 6 ^^'^ = T. From inequality 

and since <J> is upper semi-continuous, we deduce that 

Ht,y,y) < limsup<&(t"^, y?> 5X( ,y* 2 e > 5X( ) 
< $(T,y* 5 ' A >V W ), 

which implies 

vl(t,y) - V2*(t, y) ~ 2^exp (A(T - i))y 7+1 - 2C ' 



< e*(r,y*^)-t52,(r,y*^). 

Using inequality v*(T, y* s < x >£) < v 2 *(T,y* 5 ' X ^) and sending A, 5, ( — > 0+, we have 
u*(t,y) < v*(t,y), for all (t,y) G [0,T] x (0,oo). 
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-k Case 2: If the set {e > : t* e,S '^ = T} is finite, then there exists a subsequence renamed 
(£*e,<5,A,C) e such that t *e,5,XX < t. Our aim is to construct a regular function denoted V'l 
(resp. 1/^2) satisfying inequality (|4,5p (resp. (|4.6|) ). We define ipi and ^2 as follows 

+ 6 (exp (A(T - i))y 7+1 + exp (A(T - t*^))(y^y +1 ) 



+ C( 



+ 



y7 +l ( y «.*,A,C) 7 +i 



) + *(^.W >y «.W >1/ «AA,C )> (t>y) e [0)T] x (0jOo)) 



and 



-(2/1 -2/) 



- <5 (exp (A(T - t«AA,C))( y «AA,C )7 +i + exp ( A ( T _ ^7+1 



C(- 



1 



+ 



From inequalities (|6.ip and (|6.5p . we have 

< t?(* s S/) <U(y) + \K\y. 

We define V>i and ^2 as follows 



L.) _ $(t-AA,C ) y*^,A 1 C ) ^ ) A,C )) (t>y) G [0>r] x (0jOo) . 



(6.41) 



2£%) for all < y < y x A 
■01 (t, y) for all Mf' C < y < At} c 



(6.42) 



and 



Miy 



u(y) 
2 

^2{t,y) 



for all y > 2M 2 



for all < y < y A 



»2 " 2 ' 

for all Aff ,C < y < M|' C 
<5,C 



(6.43) 



M 2 y + M 2 for all y > 2M" 2 A , 



_ 1 

where y = (7|if|) 7+1 , Mi satisfies 



Mi > ^ 2-£ + |if|V v 1 



7 



dy 



(6.44) 



y 2 = (27|JC'|) 7+1 , M 2 satisfies 



dy 



2 



(6.45) 



where v 2 *rf is the right-hand derivative of V2* with respect to the variable y and M 2 = 
v*2(i ; 2M 2 '^) — 2M2M 2 ' < ' — 1. The choice of M2 is possible since V2* is convex and nonin- 
creasing ( see assumptions of the comparison theorem ) and so 

4d(** £AA ' c ,y 2 * e,a,A ' c ) < 4d(^' M ' c ,2M«) 
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and from the optimality of (t* e ' <5 ' A,< >, y™' 6 '*'*'), we have 



()u <%* d (* W, S2/r ' '")■ (6-46) 



From the definition of y 1 and inequalities (|6,4ip . (|6.44|) . we have 



v{t, y) < 2U{y) for all < y < y % 

4" 



v{t, y) < U{y) + \K\y < M x y for all y > 2M & 2 ' C 



and so one could obtain inequality (I4.5P for ipi. 

From the definition of y 2 , inequalities (|6.4ip . (|6.46p and using the convexity of V2*d) we 
have 

/ v{t, y)>^- for all < y < y_ 2 



v(t, y) > M 2 y + M' 2 for all y > 2M 



6,C 

■l an y ^_ *±v±, 

and so one could obtain inequality (|4.6p for ip 2 - 

To prove the comparison theorem, we need to derive an equivalent formulation of viscosity 

solutions an in Soner [23] Lemma 2.1. For this, we show that the control p runs along a 

compact set and the lower semi-continuous envelope of H is continuous in its arguments 

which is the object of the next lemma. We denote by p := — — . We define E' by 

J r mm diiTi J 

Ki<d 



, . 2i\4' C + l 2M 9 5 ' C + 1 ] 

P = (p t )i< t <d, 0< Pl < * s V * g V p \ (6.47) 



and the Hamiltonian H by 



H (t, y, v, ^ = M (t, y ,v,^+y(a-P + (P-J2 Pi^i)^j j , (6-48) 



Lemma 6.1 We assume that v 2 * is convex and nonincreasing. Then, we have the following 
inequalities 

min { ^r AA ' C , yT AX ' ( ) + H (t*^<, y?> S >^, v\, ^) , (6.49) 



and 



^1 /,*6,<5,A,C , «,W \ 



>o, 



mm 



< 0. 
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Proof. If there exists < Iq < d such that pi > p, then (/3 — pi5iiri)+ = 0. Let p G £ 



i=i 



be a fixed vector. If pi > -, A ? V p, < i < d , then we have 

2/1 ' ' 

Fi(ft) := ^i(t* eAA ' c ,Piyr 5 ' A,c )-^i(t* eAA ' c ,yr 5 ' A ' c ) 

= (M X - ^(t*^A,C )yr M,C ))p . y * e AA,C _ ^(^C^.**) 
Prom (|6.44p . we deduce that the function pi — > F\(pi) is non-decreasing and so 



1=1 
(i 



i=l 



:= inf Mi*^,^'^) 
pes' 

:= v^if^'^yl 6 ' 6 '^). (6.51) 

The criterion of the optimization problem (|6.51|) is continuous with respect to p, £ is 
a compact and so there exists p* solution of (|6.5ip . We consider a sequence (t k ,y k )k £ 
[0,T] x (0,oo) such that (t k ,y k ) — > (t* e > 5 - A ^, y* e,<5,A ' C ) when n goes to infinity. We denote 
by p k the optimum i.e. 

v° pt {t k ,y k ) = fx(t k ,y k ,p k ). (6.52) 

Since p k 6 £ which is compact, then along a subsequence denoted also by (p k ) k , we 
have — )• p. From the Taylor expansion formula and using the continuity of f\ in her 
arguments, we have 

h(t k ,y k ,p k ) = h(t*^,yl M ,p) + o(l) 

> vl pt (t* e ' 5 ^,yl eAX ^)+o(l) (6.53) 

From (|6.52p and (|6.53j) . we deduce that 

vf{t k ,y k ) > vf(t*^,yl M ) + o(l). (6.54) 

To obtain the converse inequality, we have 

^i pt (4,yfc) < fi(tk,Vk,Pi) 

= fi(t* e ' sx<: ,yT' sx<: ,pl) + o(i) 

= Wl p '(t* e ' 5 ' A - c ,yr AA ' c )+o(i), 
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and so lim v° pt (tk,Uk) = v° pt (t* e,S ' X ^ , yl e ' S ' X ' < '). This proves that v° pt is continuous in 



k — »oo 



H* (W^^'^i, ^j=H (W^^'^i, ^ . (6.55) 
From equality (|6.42p . the function v\—ip\ has a strict global minimum 

[0,T) x (0,oo). Using the definition of viscosity supersolutions (see inequality (|4.5p ). equa- 
tion ()6.55p and (j6.51|) . we obtain 



mm 



^(t*^,y?> 5 ^)}>0, (6.56) 



2M 

Let p G S be a fixed vector. If <0j > t ,j 2 A f V o, then we have 

y a " 

= (m 2 - ^(t«'^«, y r AA,c ))p^r' 5 ' A ' c - ur^,^) 

From the definition of ^2 and inequality (|6.45p . we deduce that the function pi — > i^O^i) 
is non-decreasing and so 

= mf I a? y r- 5 ' A,c , ^ ^ + y r AA>C ^-/3+(/3-g ftW + j I , 

From equality (l6.43p .Ci9* — tp2 has a strict global maximum 

(0,oo). Using the definition of viscosity sub-solutions (see inequality (|4.6p ) and (|6.57p . we 
obtain 

min {^(T 8 **, y?' s ' x '<) + y? AX ' c , ^ 

-^(^■W^-AA^j^o (658) 

From inequalities (|6.56p and (|6.58p . using the fact the control p runs along a compact set 
and repeating arguments of Soner [23] Lemma 2.1, we easily obtain an equivalent formula- 
tion of viscosity solutions given by inequalities (|6.49p and (|6.50p . □ 
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We come back to the proof of the comparison theorem. Remarking that min{<i, e} — 
min{/, g} > implies either d — / > or e - <? > 0, inequalities (|6.49p and ()6.50p im- 
ply 

- 5X exp (A(T - t* 6 ' 5 '^)) ((yr' 5 ' A ' C ) 7+1 + (vT^)^ 1 ) + ?i - T 2 > 0, (6.59) 



or 



-S exp (\(T - ((yT' 5 ' X ' C V + {y* 2 tAU V ) (6-60) 

+ c ( (y r' 5 - A V 2 + (yr AA V 2 ) 



where 

r d 



Tj := inf \ Y d *i(n(t^ M ,PiV?' 8,XA ) ~ m* £ ' S ' U ,y*i e ' S ' X ' C ) (6-61) 



. i=l 



( Pi - l)yl^(-(yl^ - y^) + 5( 7 + 1) exp (A(T - t*.W))( y *M,A,C)7 

C(7 + + ^ w (« " + (/5 " E M*) + ) } , 



T 2 := inf ^^(^(^^^^-^•^^'WS''^) 



8=1 



( Pl - i) y r' 5 ' A ' c (^r AA ' f - y?' 5A< ) - % + 1) exp (a(t - t* e -^-f))(yr AA,c ) 7 



+ C(7 + l) ' 



and £ is defined in (|6.47|) . From Lemma 15. II vo* is continuous w.r.t the state variable (£2* 
is convex on (0, 00)) and so the criterion of the optimization problem (|6.6ip is continuous 

w.r.t p. The set £ is compact and since lim V2*(t,y) = 00, there exists a solution denoted 

y — >o 

by p* e ' S ' X '^ 7 to the optimization problem (|6.6ip satisfying p* e,s,x,< ' > 0, for all 1 < i < d. 
We define / as follows: 

m ■■= E^(^(^' 5 ' A * c ^^r AA ' ? ) -^(r eAA ' c ,yr AA ' c ) 

i=l 

- ( ft - i)yr' 5,A>c (f(j/r AA ' c - 2/r AA,c ) - + 1) ex P (a(t - r^))^) 7 

+ C(7 + i) ( ^,J, C)7+2 )) + ^r AA ' c (* - + 08 - E MirO+J , 
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From the definition of /, we have 

d 

i=i 

< C 5 ' c , (6.62) 
where 1 denotes a M d - valued vector with all components equal to 1. Since 

for all 1 < i < d, inequality ()6.59j) implies 

- 5X (exp (A(T - t«AA,C))( y «^,A,c )7 +i + exp (A(r _ t *eM,C^( y *e,S,\Cy+A > T 2 - 7\ > T^ E ' M ' C , 

where 



TP :=\ ^^((2(^-1) + ! 



2x (j/l ~ V_2 

Pi ) 

. 1=1 



+ J(( 7 + l)(p 4 - 1) - pj +1 + 1) (exp (A(T - t«AA,C ))(yr AA,C ) 7+i 
+ exp (A(T - t «.W))( y «.W)7+i N 



+ c((7 + 1)(1-P»)-4t + i1 1 - + 



/ v>r w ) 7+1 (yr AA,c ) 7+1 , 

+ (y 2 * e '* ,A,c - 2/r* 5 ' A ' c ) {a-p + ifi-Y, PiSmh) ) } • 

i=l J 

Sending e — )• + , we obtain 

- 2,5 A exp (a(T - i* 5 ^)) ( y *Wyr+i 

> X>( 2 <* (( 7 + l)( P f X( - 1) + (pf ' A ' C ) 7+1 - l) exp (A(T - t* s > x <<j) ( y *^y+i 



i=l 



+ 2(^( 7+ l)(l-pf.«)-^^ + lj^j. (6.63) 

Using the boundedness of y* <5 > A 'f and p* S ' X ^ £ £ , along a subsequence (y* S ' X ' < ', p* 6 ^^) 
converges when A — s- oo. Let's denote (y*" 5 '^, p* 5 '^) its limit. Sending A — > oo in inequality 
(j6.63|) . we obtain 



£ *i fa ((7 + l)(pf ' C " 1) + (pf ' C ) 7+1 " l) exp (A(T - f 



i=i 



*<5)C N \7+l 
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which implies, there exists iq, 1 < i$ < d such that p*^ = 0. Sending e — > and A — > oo 
in inequality f|6.62j) . we obtain /(p* 6 ''') = oo < C" 5 '^ which is false. 
Sending e — > oo in inequality (|6.60p . we have 

-5exp (A(T - I***)) (y^C) + C (r ^ C)7+2 > 0, 



which implies 

Jexp (A(T - i* 5 ^)) 



< ( 6 - 64 ) 



C ~ (Mi ,c ) 2 (-v+i) ' 

Using the boundedness of t* 5 ' A '^, along a subsequence t* 5 ' A >£ converges when A — > oo. From 
inequality ()6.64p . we have necessarily i*" 5 ^ = T. 
From inequality 

$(t,y,y)<<I>(r^,yr' 5 ' A * C ,yr' W ) 
and since 3> is upper semi-continuous, we deduce that 

$(t,y,y) < limsuplimsup$(t* £ ' 5 ' A ' C ,y* eAA ' C ,^ e ' 5 ' A ^) 

A— >oo e— >0 

< vUT,y S ' C )-v 2 *(T,y s '<)<0. 
Sending 5, (, — > + , we obtain 

v{{t,y) < v 2 *(t,y), for all (t,y) G [0, T] x (0,oo). 
and so Theorem 15.11 is proved. □ 
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